Dephasing and Metal-Insulator Transition 



O 

o 
o 

(N 
o 

Q 



■4— » 
•4-* 



o 
o 



> 
in 

(N 
(N 

o 
o 

-I— > 
c3 



o 
o 



X 




Tl 




T2 






Junren Shi and X. C. Xie 

Department of Physics, Oklahoma State University, Stillwater, OK 74078 

The metal-insulator transition (MIT) observed in two-dimensional (2D) systems is apparently 
contradictory to the well known scaling theory of localization. By investigating the conductance of 
disordered one-dimensional systems with a finite phase coherence length, we show that by changing 
the phase coherence length or the localization length, it is possible to observe the transition from 
insulator-like behavior to metal-like behavior, and the transition is a crossover between the quantum 
and classical regimes. The resemblance between our calculated results and the experimental findings 
of 2D MIT suggests that the observed metallic phase could be the result of a finite dephasing rate. 

PACS numbers: 71.30.+h, 73.40.Hm 



Since the discovery of the metal-insulator transition 
in two-dimensional (2D) systems |J, several theoretical 
models have been proposed to understand the phenom- 
ena. Among the proposed theories, some |^-||] can be 
considered as semi-classical theory. Although they are 
different in details, the basic idea is the same. They all 
consider that the metallic phase observed in the exper- 
iments is the classical phase, i.e., the metallic behavior 
of the system can be well understood under the classi- 
cal picture. For instance, the percolation model, initially 
roposed by He and Xie ^ and further extended by Meir 
provides good description of many experimental facts. 
In this approach, the system consists of inhomogeneous 
carrier distribution with high density conducting regions 
and low density insulating regions. By assuming that the 
electrons are totally dephased in each separate region, 
one can consider the metal-insulator transition as the 
classical percolation transition of the high density con- 
ducting regions, and calculate the total conductance of 
the system by classical random resistance network. The 
dephasing of the carriers is essential for the model since 
a pure quantum system will never percolate according to 
the well-known scaling theory of localization || . 

It is not obvious that a system can be considered as 
a classical system at low temperatures. A recent experi- 
ment |jj on 2D systems shows that the phase coherence 
length is quite long, typically 600-1000 nm. Thus, it is 
more likely that a real system is in the regime where 
quantum effects compete with classical effects. Classical 
effects are manifested by a finite phase coherence length, 
which may be due to a finite temperature or other novel 
mechanisms. For instance, there are experimental indi- 
cations that the dephasing rate may be finite even at zero 
temperature ||. On the other hand, in a disordered 2D 
system, the quantum effect always causes the localization 
length to be finite. For a 2D system with metal-insulator 
transition, the localization length strongly depends on 
the carrier density. Actually, by changing the carrier den- 
sity, the conductance of the system may change several 
orders JlJ which implies a substantial change in the local- 
ization length. The behavior of the system is determined 
by two competing length scales: localization length and 
phase coherence length. Therefore, it is important to 



FIG. 1. The model to simulate the system with dephasing 
and disorder. The rectangles represent the normal random 
scatters, and the diamonds represent the dephasing scatters. 



study transport properties by varying these two lengths. 

In this paper, we study the interplay of disorder and 
dephasing of one-dimensional (ID) systems in transport 
properties. We limit to ID models to reduce the se- 
vere finite-size effect in numerical results for higher di- 
mensions, although our conclusion can be carried over 
to higher dimensions. Our ID model (Fig.l) consists of 
normal random scatters and dephasing scatters, alterna- 
tively. While the normal random scatters give rise to a 
finite localization length, the dephasing scatters random- 
ize electron phase. The normal random scatters are con- 
structed by M 5-barriers with random height qi, which 
has the distribution 
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In the model Aq controls the randomness of the system. 
The transmission and reflection coefficients for the nor- 
mal random scatters can be calculated from the transfer- 
matrix for individual i5-barrier 101, 
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where t and r are the transmission and reflection am- 
plitudes for the barrier, and k is the momentum of the 
injected electron. The transfer matrix for M sequential 
(5-barriers can be calculated from 
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where X is the transfer matrix describing the propaga- 
tion of the electron from one 5-barrier to the next. As- 
suming the spacing between the neighboring barriers is 
unity, X is 
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FIG. 2. The structure of an individual dephasing scatter. 
The system is connected to two identical external electron 
reservoirs PI and P2. An electron can be scattered along the 
directions shown by the arrows. There is no backscattering 
caused by this dephasing scatter. 
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Using the transfer-matrix technique, the localization 
length can be determined analytically || . 

To introduce the dephasing effect into the system, one 
has to include the interaction between the system and the 
environment. The Biittiker model [ fio[ shown in Fig.^ is 
the simplest way to achieve that. In this approach, the 
system is connected to the external electron reservoirs 
via the dephasing scatters. With a certain possibility, an 
electron is scattered into the external reservoirs, totally 
losing its phase memory, and then re-injected into the 
system. Two restrictions are imposed to reflect physical 
reality. First, the net current between the system and 
the reservoirs should be zero so that each scattered elec- 
tron will finally return to the system. To do so, one can 
adjust the chemical potential fj, of the external reservoirs 
such that I pi + Ip2 = 0, where Ip\{Ip2) is the current 
between the system and the external reservoir PI (P2). 
Second, the system is connected to two identical electron 
reservoirs PI and P2, and the S- matrix between the sys- 
tem and the reservoirs is designed so that the electron is 
only scattered forward, thus the dephasing scatters will 
not cause any momentum relaxation. The 5-matrix reads 
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where a is the possibility that the electron is scattered 
to the reservoirs, namely, the dephasing rate. The phase 
coherence length is estimated by L v ~ M/a. 

The localization length of the system is determined by 
the normal random scatters, while the phase coherence 
length is determined by the dephasing scatters. For a sys- 
tem with N dephasing scatters, there are N + 2 external 



chemical potentials fii with i = L, R, 1, 2, • • • N, where 
fJ-L(p) is the chemical potential for the left (right) mea- 
surement electrode. The system satisfies the multi-lead 
Ohm's law, 

h = ViiMi h 3 = L, R, 1, 2, • • ■ N. 



The conductance between the leads, cr.y, is determined 
by the Landaur-Biittiker formula, 
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where Xy is the transmission coefficient between the 
leads. For the leads attached to the dephasing scatters, 
the transmission coefficient is the sum of all transmission 
coefficients between leads PI and P2, 

Tij = Ti pi j pi + Ti pi j P2 + Ti p2j pi + Ti P2j P2 , 
foH,j = l, 2, ---N. 

The gauge invariance, namely shifting each /i, by a con- 
stant should not affect the result, is satisfied through the 
condition 

era = -22<Tij. 

The total conductance of the system is calculated by im- 
posing the condition 

l L = fori = 1, 2, ■ ■ ■ N. 

After some algebra, the total conductance of the system 
can be written as 
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Although we start from a seemly artificial model, the 
total conductance formula actually reflects physical real- 
ity. It contains two kinds of contributions: the first is 
the direct conductance ulr coming from direct quantum 
tunneling; the second is the correction due to the dephas- 
ing effect caused by classical sequential tunneling. Thus, 
the conductance formula is consistent with the general 
picture of the dephasing effect on a conductance. 

The typical behavior of the conductance for the sys- 
tem is shown in Fig.^, where we plot the conductance as 
a function of the dephasing rate a. In a real system, the 
dephasing rate is a monotonic function of temperature, 
so the plot can also be considered as the temperature 
dependence of the conductance. We have systematically 
calculated the conductance for different sets of parame- 
ters, and the results show qualitatively similar behavior, 
although the peak position depends on the parameters. 
The most important feature of the plot is that the con- 
ductance is not a monotonic function of the dephasing 
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FIG. 3. The typical behavior of the conductance. We use 
the parameters Q = q/2k = 1.1, Aq = 0.5, M = 5, N = 100, 
and average over 500 samples. The gray and white regions 
show the quantum and classical regime respectively. 
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FIG. 4. The comparison between the localization length 
and the l/a pea k, which is equivalent to the phase coherence 
length. Qo = qi/2k, and Aq = 0.5. 

rate, or temperature. The plot can be divided into two 
regions, the gray and white regions. In the gray region, 
the dephasing rate is low and the phase coherence length 
is long. The electron is localized within the phase coher- 
ence length, so the system shows quantum localization 
with insulator-like behavior, namely, the conductance de- 
creases as temperature drops. This region can be consid- 
ered as the quantum region. On the other hand, when 
the dephasing rate becomes higher, the system enters into 
the classical regime and the conductance increases with 
decreasing temperature, a typical metallic behavior. 

The turning point between the quantum and classi- 
cal regimes can be determined by comparing the phase 
coherence length with the localization length £ calcu- 
lated from the transfer matrix formalism The result 
is shown in Fig.^. The phase coherence length L v is ob- 
tained by using the value of a corresponding to the peak 
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FIG. 5. The dephasing rate (or temperature) dependence 
of the resistance for different carrier densities (Qo = qi/2k). 
The gray region shows inaccessible region due to the finite 
cutoff of the dephasing rate. We use Aq = 0.5, and Qo =0.1, 
0.5, 1.0, 1.1, 1.2, 1.3, 1.5, 2.0. 

in the conductance plot (see Fig. 3) for a given Q . One 
can see that the two lengths are approximately equal. 
Thus, it shows that the transition occurs at the point 
where L v ~ £. From Fig.|| and Fig.^J, one can clearly 
see that by changing the phase coherence length or the 
localization length, it is possible to observe the transi- 
tion from the metal-like behavior to the insulator-like 
behavior, and the behavior of the system is determined 
by the competition between the localization length and 
the phase coherence length. When L v ~ £, the quantum 
physics of localizaiton ceases to exist. Therefore, we be- 
lieve that it is the transition between the quantum and 
the classical phases. 

Great resemblance can be found when comparing the 
experimental findings m with the results from this sim- 
plified model. In Fig||we plot the dephasing rate (or 
temperature) dependence of the resistance. Different 
curves are for different scattering parameter Q = q/2k. 
Changing the carrier density is equivalent to changing 
the Fermi momentum, which in turn changes Q. Thus, 
different curves correspond to different carrier densities. 
In the plot, we impose a finite cutoff of the dephasing 
rate a, which makes lower a inaccessible (gray region in 
the plot). For a above the cutoff, the phase coherence 
length is always finite. Depending on whether is 
larger or smaller than the localization length £, one can 
either observe the metal-like or the insulator-like behav- 
ior, as shown in FigJ|. When the finite cutoff falls upon 
the turning point (the maximum conductance point in 
Fig. 3), the system shows the "critical" behavior, where 
the resistance is nearly flat within a certain temperature 
range. In Fig|| we show the density dependence of the 
conductance for different a (or temperature). One can 
identify a "fixed" point where different curves cross at. 
Similar feature has been seen in many experimental plots. 
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FIG. 6. The density dependence of the conductance for 
different temperatures. We reverse the direction of the axes 
so that it is easier to compare with n — p plot wildly used in 
the literatures. 



The finite cutoff of the dephasing rate can be justified 
by two possible reasons. First, the cutoff may be due to 
a finite temperature. In this case, if we assume the de- 
phasing rate goes to zero when temperature approaches 
zero, as suggested by a simple power law a ~ T u , there 
is always an upturn of the resistance at low enough tem- 
perature as shown in Fig. 5 for low values of a. This sug- 
gests a re-entrance to an insulator at low temperatures. 
The re-entry behavior may have already been observed 
in the recent experiment Under the circumstance, 

the transition is a finite temperature effect. The second 
possibility is that the dephasing rate might be finite at 
zero temperature ||. Consequently, the metallic phase 
will survive even at zero temperature. If we adhere to 
the original definition that a metal has a finite resistance 
at zero temperature while the resistance of an insulator 
diverges, the system will always be a metal. The reason 
is that on the low density "insulator" side, the resistance 
will saturate to a finite value at T = because of the 
finite phase coherence length. However, in a similar plot 
as shown in Fig. 6, a "fixed" point can still be identified 
which can be used as an operational definition of " metal- 
insulator transition" . 

The saturation of the dephasing rate at low tempera- 
tures is still a controversial issue. Some argue that the 
saturation observed in the experiments is not an intrinsic 
effect. Nevertheless, whether it is intrinsic or extrinsic, 
the same factors which cause the saturation should have a 
similar effect on the conductance. To justify a real metal- 
insulator transition, one has to clearly rule out those ex- 
ternal factors that may cause finite dephasing rate at low 
temperatures 

In summary, we have studied the interplay between 
dephasing and disorder. Based on a ID model, we show 
that by changing the phase coherence length or the lo- 
calization length, it is possible to observe the transition 



from the insulator-like behavior to the metal-like behav- 
ior, which corresponds to a transition between quantum 
and classical phases. The great resemblance between the 
results from this simplified model and the experiments 
suggests that the quantum effect is important at low 
temperatures, although the high temperature behavior 
is dominated by the classical effect. We suggest that 
conductance experiment should be accompanied by a de- 
phasing rate measurement to address the effect of a finite 
coherence length. Although our calculation is for a ID 
model, the same physics should survive at higher dimen- 
sions. 
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